Magnetic and pair correlations of the Hubbard model with next-nearest-neighbor 

hopping. 



Alain F. Veilleux, Anne-Marie Dare, Liang Chen, Y.M.Vilk, and A. -M.S. Tremblay 
Departement de physique and Centre de recherche en physique du solide. 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1 



(13 July 1995, sond-mat/950704! 



A combination of analytical approaches and quantum Monte Carlo simulations is used to study 
both magnetic and pairing correlations for a version of the Hubbard model that includes second- 
neighbor hopping t' = — 0.35t as a model for high-temperature superconductors. Magnetic properties 
are analyzed using the Two-Particle Self-Consistent approach. The maximum in magnetic suscep- 
tibility as a function of doping appears both at finite t' and at t' = but for two totally different 
j^-^ | physical reasons. When t' = 0, it is induced by antiferromagnetic correlations while at t' = — 0.35t 

■ it is a band structure effect amplified by interactions. Finally, pairing fluctuations are compared 

^\ ' with T-matrix results to disentangle the effects of van Hove singularity and of nesting on supercon- 

7—i , ducting correlations. The addition of antiferromagnetic fluctuations increases slightly the d-wave 

• superconducting correlations despite the presence of a van Hove singularity which tends to decrease 

them in the repulsive model. Some aspects of the phase diagram and some subtleties of finite-size 
scaling in Monte Carlo simulations, such as inverted finite-size dependence, are also discussed. 

^t" ■ 

t-H ' PACS numbers: 75.30.Cr, 74.72.-h, 75.10.Lp, 75.40.Mg 

i— I ! 
> ■ 

^ ' I. INTRODUCTION 

£ — , High-temperature superconductors pose a theoretical challenge pj for many reasons. One of them is that their band 
structure does not suffice by itself to explain most experiments, even in the normal state. This is one of the reasons 
that much work has been devoted to the simplest model which incorporates the effects of short-range interactions, 
the Hubbard model, 

-i— » ; 

SH = - ]T t itj (c\ a c ja + c] a a^j + U Y^ n it n H ■ (!) 
<ij>(T i 

In this expression, the operator ci a destroys an electron of spin a at site i. Its adjoint c] CT creates an electron. 
The symmetric hopping matrix £jj determines the band structure. We consider the case where tjj = t for nearest- 
O ■ neighbors, ti^ = t' for next-nearest neighbors and £y = otherwise. Double occupation of a site costs an energy U 
due to the screened Coulomb interaction. We work in units where the lattice spacing is unity, ks = 1, h = 1 and 
t = 1. Since high-temperature superconductivity generally occurs in systems with a pronounced planar structure, the 
two-dimensional square lattice version of the model is considered, as usual. 

The nearest-neighbor (t 1 = 0) version of the above Hubbard model has been extensively studied. Within this model, 
the only hope to get even qualitative agreement with experiment is if the interactions or the antiferromagnetic fluctu- 
ations are very strong. Indeed, with the nearest-neighbor model the Fermi surface topology and the filling dependence 
of both the Hall coefficient and the uniform magnetic susceptibility are all qualitatively wrong, a disagreement that 
cannot be removed perturbatively. On the other hand with second nearest-neighbor hopping, the band structure 
becomes more realistic and all the above physical quantities as well as the position of neutron scattering intensity 
maxima [|| (D @] Hi have at least the correct qualitative behavior. 

It is thus important to understand the Hubbard model Eq.(Q) including both nearest-neighbor hopping t and 
next-nearest-neighbor hopping t' since weaker interactions in this more realistic model might suffice to explain the 
experimental data. Two main questions are addressed here by a combination of numerical @ @ |§| [|| JTo| 0] JlJ] 
and analytical approaches: jb| First, the origin of the maximum in the uniform magnetic susceptibility as a function 
of band filling at low temperature, and second the relative influence of density of state effects (van Hove singularity) 
and Fermi surface-topology effects (nesting) on pairing correlations. While both van Hove singularity and nesting 
occur simultaneously at half-filling in the usual t' = model, there is no nesting when t' ^ 0. |14| 

In the rest of this paper we first discuss the Monte Carlo approach and the range of parameters studied. In this 
same section, we discuss subtleties of finite-size scaling including a new effect, inverted finite-size scaling. In the 
following two sections, we present results in turn for the magnetic susceptibility and for antiferromagnetic as well as 
pairing correlations along with analytical |15| approaches to understand the results. It will be shown that for magnetic 
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properties, the recently proposed Two-Particle Self-Consistent (TPSC) approach |L6| compares very well with Monte 
Carlo results, allowing us to understand the origin of the maximum in magnetic uniform magnetic susceptibility as 
well as some features of the phase diagram. The same kind of agreement was found before in the case of the nearest- 
neighbor model. Jl6| Superconducting fluctuations on the other hand are accounted for by T-matrix effects far from 
the van Hove singularity |L7j E8| . The physical interpretation of the results is summarized in the conclusion. 



II. QUANTUM MONTE CARLO APPROACH, SIMULATION PARAMETERS AND FINITE-SIZE 

EFFECTS 

After a brief discussion of technical details of the simulations, we make several general points about finite-size 
effects. These effects can be especially unusual in the presence of t' . We point out in particular how false signals of 
antiferromagnetism can be detected as an inverse finite-size dependence. Finite-size effects will be discussed further 
in the other sections for each specific physical situation encountered. 



A. Simulations 



We use the so-called determinantal quantum Monte Carlo approach |2(J |21| |2^] p3| . As usual the sign of the 
fermion determinant renders the statistical accuracy very poor in certain regions of parameter space so that long com- 
putation times are required. Figure 1, computed using the standard gauge 24 for the Hirsch-Hubbard-Stratonovich 
transformation, gives the average sign of the fermion determinant for typical parameters considered in the present 
paper. For t' = 0, the sign was always positive at half-filling because of particle-hole symmetry [|25| but this is no 
longer the case for t' ^ 0. 

All the simulations presented below are for either t' — 0, or t' = —0.35 which, in addition to minimizing finite-size 
effects, as described below, is typical for high-temperature superconductors. Indeed, band-structure calculations |2(| 
|3 suggest t' = -0.16 for La 2 CuO A and t' = -0.45 for YBa 2 Cu 3 7 . 

Two to five different runs were done for each point with 6 x 10 5 to 2 x 10 6 measurements per run. Measurements 
were grouped into 400 equal-size blocks to estimate the statistical error. Between 10 3 and 10 4 flips per spin in space- 
time were performed for the warm-up period. Sparse matrix techniques were used for the exponential of the kinetic 
energy operator. |^3| The most difficult points took 300 hours each to be obtained at a computing speed of 26 million 
floating point operations per second. The imaginary-time discretization step was between 1/8 and 1/16. Some results 
were also obtained on a Fujitsu supercomputer. 



B. Finite-size effects 



One should recall that finite-size corrections arise at low temperature because of two different physical effects: a) 
The correlation length arising from interaction-induced collective phenomena (antiferromagnetism) can become larger 
than the system size, b) The temperature can become smaller than the energy-level separation or, equivalently, the 
thermal de Broglie wavelength £ f /j = vpj (ttT) can become larger than the system size. In the latter case, even the 
results for U = contain finite-size effects. 

Since finite-size effects limit the lowest temperature that we can consider, they also determine the smallest accessible 
value of the next-nearest neighbor hopping t' . Indeed, non-trivial effects caused by t' appear only when t' > T. 
Furthermore, it is important to realize that finite-size corrections may depend on filling at fixed temperature. This 
is particularly important in the context of one of our problems, namely the maximum in the uniform magnetic 
susceptibility as a function of doping at low temperature. Our results, in Figs. 2 and 3, will be discussed in more detail 
in the next section. To understand finite-size effects here, we first recall that the uniform magnetic susceptibility of 
the non-interacting infinite system has a maximum away from half-filling as long as t' ^ 0. In a 10 x 10 system at 
U = 0, P = 1/T = 5, one can check that this maximum first appears away from half-filling only for t' in the vicinity 
of —0.35. Smaller values of t' do not even lead to a filling dependence qualitatively similar to that of the infinite 
system at zero-temperature. For t' = —0.35, several of the allowed wave vectors of the finite system are close to the 
Fermi surface at the filling n about 0.6 where is located the van Hove singularity in the infinite-size system. This 
leads to a better finite-size representation of this singularity. The finite-size effects in the non-interacting 8x8 system 
illustrated by the dashed line in Fig. 3 lead to the small oscillations around the infinite-size result. These oscillations 
would become much larger at lower temperature where they could lead to erroneous conclusions. 
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The interplay between finite-size effects and antiferromagnetic fluctuations in the finite t' model is also very peculiar. 
Antifcrromagnetic fluctuations may be observed in the (tt, tt) component of the magnetic structure factor defined by, 

s ( q ) = ]v ^ exp ^ ~ r ^ ^ ni,T ~ ^ njA ~ nu ^ ' 

Fig. 4 allows one to suggest that in the t' = —0.35 model the antiferromagnetic fluctuations are largest around n = 1 
since this is where finite-size effects are largest. By contrast with the t' = model however, we observe here an inverted 
finite-size effect. What we mean is that the magnetic structure factor is converging down to the infinite volume result. 
In the t' = model, this is never observed. When the correlation length is finite but larger than the system size at this 



temperature, S (q) normally increases with system size until it saturates. The TPSC described below |16 reproduces 
the inverted finite-size effect of the t' = —0.35 model and suggests that Monte Carlo calculations on 10 x 10 would 
be necessary to converge to the infinite-size result at this temperature. The source of the inverted finite-size effect 
can be found already from the U = calculation which shows that for t' = —0.35 the magnitude of the q = (tt, tt) 
susceptibility ona6x6ora8x8 system is anomalously large compared with the infinite-size limit. This can lead to 
spurious results in zero-temperature calculations. Indeed, in Ref. |1I[ it was found that the 6 x 6, t' = ±0.4 system at 
half-filling has long-range order at zero-temperature. Our analytical TPSC approach suggests that long-range order 
is present at T — on both 6x6 and 8x8 system but not in the infinite volume limit so that the disappearance of 
long-range order at zero temperature would be apparent, in this case, only for very large systems. 



C. Finite-size analytical calculations 

The above considerations show that when an analytical approach |l5[ is available, it can be used to understand 
finite-size effects by simply comparing the infinite-size limit with calculations on systems that have the same size as 
those where the simulations are done. It should also be clear that even when an analytical approach is not available, the 
non-interacting case can carry useful information on spurious finite-size effects that arise when the thermal de Broglie 
wave length is larger than the system size. A striking example of this is provided by Fig. 5. At finite temperature 
((3 — 5), the non-interacting susceptibility has qualitatively the same overall magnitude and wave vector dependence 
on a very large lattice (Fig. 5a), and on a small 8x8 lattice (Fig. 5b). In the zero-temperature limit however, Figs. 
5c and 5d show that the 8x8 lattice is qualitatively different from the infinite-lattice result. In particular, on a finite 
lattice the maximum susceptibility is at zero wave vector and is about five times larger than the true maximum that 
is furthermore located at a totally different wave vector. 

Analytical calculations on finite lattices also have difficulties of their own in the presence of interactions. For 
example, the Random Phase Approximation (RPA) can be more divergent on a finite lattice at finite temperature 
than the corresponding calculation on an infinite lattice. Indeed, consider the Lindhard function 

Xo(q,^ = 0) = 4E /(£k) ' /(£k+q) ( 2 ) 
k Ck _ £k+q 

where the sum is over the N discrete lattice points of the Brillouin zone, and / (ek) is the Fermi function. Suppose 
that for to of the wave vectors on the Fermi surface, the condition ek = fik+q is satisfied. Then the total contribution 
of these points to the Lindhard function is given by 



to <9/(e k ) 



N 8e k 



ATN 



(3) 



On a finite lattice this contribution becomes huge at small temperatures T <C to/ (AN) while it is never the case in the 
infinite-size N — > oo limit. Hence, the RPA susceptibility xo/ 0- ~ Uxo) may be divergent on a finite lattice at low 
temperature, even when for an infinite system the condition Uxo I ensuring the validity of RPA is satisfied. Since 
susceptibilities are always finite on finite systems, this shows that perturbative techniques that are appropriate for 
infinite lattices may fail when applied on small systems. In the case just discussed, the infinite-lattice result obtained 
from RPA may be closer to the true finite-lattice result than the same RPA calculation done on the finite system. 
Consequently, in the TPSC approach |Q described below, the intermediate quantity xo (q, iln = 0) appearing in the 
calculations sometimes has to be computed for a discrete set of q wave vectors but with a sum over k done for the 
infinite system to correctly estimate finite-size effects. Note that when the de Broglie wave length is smaller than the 
system size, then the inequality (£th = «f/ (kT)) < N 1 / 2 ensures that the troublesome condition T <C to/ (AN) does 
not occur since T > vf/ (ttN 1 / 2 ) generally implies that T 3> 1/N is satisfied. 
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III. MAGNETIC STRUCTURE FACTOR AND TWO-PARTICLE SELF-CONSISTENT 

APPROXIMATION. 



The Two-Particle Self-Consistent approximation (TPSC) [pl| agrees at the few percent level with simulations on 
the nearest- neighbor model up to intermediate coupling (U < 8). The agreement is overall better than FLEX or 
parquet. [^8| While direct diagrammatic approaches [^9| |30| can work well far from half-filling, they fail as soon as 
antiferromagnetic fluctuations start to increase. Hence, we compare our Monte Carlo results only with the TPSC 
approach. We first briefly discuss the approach and use it in following subsections to understand the Monte Carlo 
results. 



A. TPSC approach 



The TPSC approach |16| |31[| can be summarized as follows. One approximates spin and charge susceptibilities 
Xspj Xch by RPA-like forms but with two different effective interactions U sp and U c h which are then determined 
self-consistently. Although the susceptibilities have an RPA functional form, the physical properties of the theory 
are very different from RPA because of the self-consistency conditions on U sp and U c h- The necessity to have two 
different effective interactions for spin and for charge is dictated by the Pauli exclusion principle (n 2 a ) = (n a ) which 
implies that both Xsp and Xch are related to only one local pair correlation function {n^n{). Indeed, using the 
fluctuation-dissipation theorem in Matsubara formalism and the Pauli principle one can write: 

— y Xch(q) =n + 2(n T n,) - n 2 = — V ^ , (4) 



,(q) = n- 2(n T nx) 



—y- 

BN ^ 1 



xo(q) 



<Xo(<7) 



(5) 



where B = 1/T, n = (n-\ ) + («x), q = (q, iq n ) with q the wave vectors of an N site lattice, iq n Matsubara frequencies 
and xo(q) the susceptibility for non- interacting electrons. The first equalities in each of the above equations is an exact 
sum-rule, while the last equalities define the TPSC approximation for Xch(q) and for Xsp{q)- In this approach, the value 
of (n^ni) may be obtained self-consistently jl(| b y adding to the above set of equations the relation U sp — ^ with 

<7H (0) = (ri|ri|) / (rij) (ri|) . As shown in Ref. |p_6| , the above procedure reproduces both Kanamori-Brueckner screening 
as well as the effect of Mermin- Wagner thermal fluctuations, giving a phase transition only at zero-temperature in 
two dimensions. In general, there is however a crossover temperature Tx below which the magnetic correlation length 
£ grows exponentially. Quantitative agreement with Monte Carlo simulations is obtained Jl£] for all fillings and 
temperatures in the weak to intermediate coupling regime U < 8. The equation for charge Eq^p) is not necessary to 
obtain the spin structure factor discussed below. 



B. Numerical results and interpretation of the maximum in the uniform magnetic susceptibility 

As above, at q = the magnetic structure factor is related by the fluctuation-dissipation theorem to the uniform 
magnetic susceptibility, x = Xsp (q — 0) which can be measured either directly |32| or through NMR Knight shift. 
One of the puzzling aspects of this quantity in high-temperature superconductors is the doping dependence of its 
extrapolated zero-temperature value. It increases with doping before decreasing again for larger dopings. 

The uniform magnetic susceptibility x that we obtain numerically is plotted as a function of filling in Fig. 2 for 
the nearest-neighbor model (hopping parameter t) J33| p4| and in Fig. 3 for the model where next-nearest neighbor 
hopping (hopping parameter t') is allowed. Recalling that filling (n), that measures the average number of electron per 
site, is equal to one minus doping, the experimentally observed maximum is qualitatively described by both models. 
The physical origin of both results is however very different. In the t' = model, Fig. 2, the maximum comes purely 
from the interactions since the U = model does not have a maximum as a function of doping. By contrast, the 
model t' =/= in Fig. 3 has a maximum even at U — (dashed-dotted line). This maximum is simply enhanced by 
interactions, as we now discuss. 

In the t' case, Fig. 3, the magnitude of finite-size effects coming from the thermal de Broglie wavelength can 
be estimated from the {7 = result that is plotted for both 8x8 and infinite lattice. Figs. 6a and 6b show that at all 
q the overall agreement is excellent between the TPSC approach iQ and Monte Carlo simulations of the structure 
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factor. The two fillings illustrated correspond to half-filling and to a filling close to the van Hove singularity. We have 
already commented on the discrepancy for the half-filled case at q = (tt, tt) in the section on finite-size effects. For 
other q points, the discrepancy with the TPSC is larger for the q = results as in Fig. 3 because in TPSC it is only 
the integral over all q that is determined self-consistently. Given the difference in scale between the results at q = 
and q = (tt, tt), it is understandable that the q = result is less precise. Finite-size analysis of the U — case and 
TPSC itself (finite mesh calculation) suggests that Monte Carlo calculations on 10 x 10 would suffice to completely 
eliminate finite-size effects. The TPSC then allows us to conclude that the physical origin of the maximum in \ 
as a function of filling for the t' = —0.35 model in Fig. 3 is a band-structure effect (van Hove singularity) which is 
already present for the non-interacting system (U = 0) and only enhanced by interactions. No complication in the 
interpretation can arise from antiferromagnetism since, as discussed in the section on finite-size effects, it does not 
occur in this t' — —0.35 model even at zero temperature. 

By contrast, in the t' = model of Fig. 2 the mechanism by which interactions create the maximum is not so clear. 
A phenomenological interpretation in weak coupling has been proposed. |35| ] However, the TPSC approach does not 
show a maximum nor does the slave-boson approach for this value U = 4. This qualitative disagreement strongly 
suggests that the maximum is a self-energy effect (appearance of a pseudo- gap ). This self-energy effect is important 
in strong coupling where both slave-bosons |36| ] and numerical simulations p3[ | do show a maximum. It can also be 
important in the presence of strong antiferromagnetic fluctuations even at relatively weak coupling U — 4 but it is 
not included in the present version of either the TPSC approach or slave-bosons. However, one can show that the 
self-energy obtained at the next level of approximation in a manner consistent with the TPSC approach does show a 
pseudogap. |3l[] This pseudogap appears as soon as the antiferromagnetic correlation length becomes larger than the 
thermal de Broglie wavelength £ t /j. The pseudogap is more pronounced at half- filling, leading to a relative decrease 
of the density of states, and correspondingly of the magnetic susceptibility. This is what might explain the maximum 
away from half-filling: As one increases the electronic density towards half-filling, the magnetic susceptibility increases 
because of the increase in density of states until the presence of antiferromagnetic fluctuations creates a pseudogap 
that depletes the density of state more and more as we move towards half-filling. 

Finite-size effects on the pseudogap are negligible as long as Cth is smaller than the system size, even when the 
antiferromagnetic correlation length is larger than system size. [ pl[ The TPSC as well as numerical simulations at 
q = (tt, tt) do show that, for the parameters of Fig. 2, we are precisely in this situation, namely the antiferromagnetic 
correlation length is larger than ^ t h but ^ t h is less than or of the order of the system size. p§| ] Hence the apparent 
convergence of the Monte Carlo results as a function of system size, despite the large antiferromagnetic correlation 
length, would be consistent with the above self-energy (pseudogap) interpretation of the maximum. 



C. Antiferromagnetic correlations vs zero-temperature phase diagram 

Recently, Duffy and Moreo (l2| have found antiferromagnetism for the t' ~ —0.2 model at half-filling for both U — 4, 
(3 = 6, and U = 6, /3 = 4. They also found that antiferromagnetism disappears away from half-filling. This seems in 
contradiction with the zero-temperature phase diagram |l4[ ] which predicts that antiferromagnetism would be favored 
near the filling corresponding to the van Hove singularity that is here at finite hole doping. We checked that the TPSC 
approach accurately reproduces the results of Duffy and Moreo |l^| and allows us to understand their origin. The 
crossover to a regime where in two dimensions the antiferromagnetic correlation length starts to grow exponentially 
p6| occurs at a temperature Tx that is comparable with t' , hence the details in the band structure introduced by t' 
are not crucial and antiferromagnetism appears first at n = 1, as in the t' = case. Note that even in this t' = 



case, the TPSC predicts 16 that slightly away from half-filling the crossover to a rapidly growing correlation length 
occurs at the antiferromagnetic wave vector, not at an incommensurate wave vector, again in apparent contradiction 
with the zero-temperature phase diagram. The key point in all this is that the shape and position of the maximum in 
Xs P (<iMn = 0) at T = Tx can be qualitatively different from that at T = 0. This shows that it is a misconception to 
believe that the zero-temperature phase diagram determines the high-temperature finite-size calculations, or even the 
real three-dimensional finite-temperature phase transition. In fact whether the finite-temperature phase transition is 
commensurate or incommensurate may be predicted in a more realistic manner from TPSC or finite-size scaling at 
finite temperature rather than from zero-temperature calculations. For example, two-dimensional systems do not show 
a real phase transition until zero temperature. Nevertheless, if finite-size calculations reveal that the antiferromagnetic 
correlation length starts to grow with system size at Tx, then in real systems small three-dimensional effects would 
generally lead to long-range order at the wave vector at which the two-dimensional \sp has a maximum. Jl8| [l9| ] 
In the t' = case slightly away from half-filling this would be at the antiferromagnetic wave vector, not at the 
incommensurate one. 
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IV. SUPERCONDUCTING CORRELATIONS AND T-MATRIX EFFECTS 



The order parameter for superconductivity is denned by, 

i i,v 

where the sum over v is over the nearest- neighbor sites of i. The form factor determines the spatial symmetry (s, 
extended S,p, d) of the Cooper pair. Q For singlet order parameters the form factor is even in space g a (v) = g a (—v), 
and the spatial symmetries usually considered include: 

s-wave, with g a (v) — 5 V fi (7) 
extended S -wave, with g a (y) = 5„ iX + <5„,_ x + 8 v<y + 6 v ,- y (8) 

d-wave, with g a (y) = 5 UjX + 5, A _ X — 5 U) y — Su,-y (9) 

Our results for the next-nearest-neighbor model t' — —0.35 are shown in Figs. 8 and 9. Size effects are still quite 
large. At this temperature, we cannot really expect results which would converge to the infinite size limit before 
10 x 10 lattices. To account for this, the calculations for the T-matrix are done for systems of the same size as the 
Monte-Carlo simulations, namely 8 x 8 or 6 x 6. 

For the nearest-neighbor model (t' — 0), we have shown in an earlier paper [|l7j that most of the correlations at 
low density can be accounted for by two-body particle-particle scattering, namely T— matrix effects. Near half-filling, 
one observes deviations from T— matrix predictions, especially for the d-wave correlations. Fig. 7, reproduced from 
Ref. jTrfl , shows that the T-matrix underestimates superconducting fluctuations at half-filling. Since for d-wave, the 
T-matrix gives in fact the same result as the non-interacting case j^| Jl8| this means that the correlations measured 
by Monte Carlo are enhanced compared with the non-interacting case. This could indicate that either the van Hove 
singularity or the nesting-induced antiferromagnetism increase d-wave correlations. 

The t' = —0.35 Monte Carlo results in Figs. 8 and 9 can help disentangle the influence of van Hove singularity 
and nesting. For these figures, we learn that: a) As expected in general and as found in particular in the t' = 
case, the T— matrix approximation works best at low electron or hole filling, i.e. near (n) = and (n) — 2. b) 
The disagreement between simulations and T— matrix is largest on the (n) < 1 side in the vicinity of the van Hove 
singularity, c) The T— matrix everywhere overestimates the superconducting correlations. Even though this should in 
no way be considered a proof, this last point suggests that the finite t' model is even less favorable to superconductivity 
than the nearest-neighbor model. Earlier studies, B Q jl(| are equally divided on this point. Here, no sign of 
superconductivity is found in the size dependence of correlations functions. A plot of the spatial dependence H pi 
would however be a more sensitive test. 

As mentioned in the context of the magnetic structure factor, the parameters studied here do not correspond to a 
regime where antiferromagnetism has started to set in at half-filling, contrary to the studies of the t' = model. For 
the t' = model the only case where the T— matrix clearly underestimates the superconducting correlations is for 
d-wave pairing near half-filling fTij| p8| at a temperature where nesting and antiferromagnetic correlations are strong. 
Since the van Hove singularity located at the same filling is expected to decrease pairing correlations, as in remark 
b) above, this seems a clear indication that antiferromagnetism increases superconducting fluctuations. We have not 
found however yet that the increase is sufficient to lead to a superconducting phase transition. To confirm the positive 
influence of antiferromagnetism on superconductivity, it would be interesting in the t' ^ case to either increase U 
or decrease temperature to reach a regime where antiferromagnetism starts to show up in the model. One should 
beware however that antiferromagnetic correlations can have a positive influence on superconducting fluctuations due 
simply to a short-range effect, as argued in Ref. ]4jJ ]. 

V. CONCLUSION 

We have studied the next-nearest-neighbor Hubbard model for t' — —0.35 and t' = using both Monte Carlo 
simulations and analytical approaches. We have also discussed at length finite-size effects. The TPSC approach |l|] 
allowed us to understand the magnetic fluctuation properties, including the occurence of antiferromagnetism observed 
in previous calculations at half-filling. |12[ 
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The uniform magnetic susceptibility obtained with the t' — —0.35 model has the same qualitative doping dependence 
as the experimental results on high-temperature superconductors. Since we have found that the TPSC approach gives 
a satisfactory explanation of the spin-spin correlation functions computed by Monte Carlo simulations, we were able 
to identify the maximum in the magnetic susceptibility as mostly a band-structure effect enhanced by interactions. 
We contrasted this with the t' = model where a maximum also shows up but for apparently very different reasons 
related probably to self-energy (pseudogap) effects which can occur either at strong coupling or in the presence of 
large antiferromagnetic correlation lengths. 

For the interaction and temperature ranges studied here, superconductivity is less favored for t' = —0.35 than 
for the t' = model since the Monte Carlo results are consistently smaller than trivial T-matrix effects. However, 
antiferromagnetic correlations do not occur in this parameter range. It is interesting to note that despite the fact 
the T-matrix overestimates superconducting correlations in the presence of a van Hove singularity, in the t' = case 
the presence of nesting and concomitant large antiferromagnetic correlations at half-filling suffice to overcome the 
detrimental effect of the van Hove singularity to increase the superconducting d-wave fluctuations above the T-matrix 
prediction. This clearly suggests that antiferromagnetic correlations can have a positive influence on d-wave pairing. 
However, the effect observed to date is clearly a short-range one so that the issue of long-range superconducting order 
caused by antiferromagnetic correlations is not settled yet. 
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Figure 1: Average sign of the determinant as a function of filling for the next-nearest-neighbor hopping model. 

Figure 2: Maximum in the uniform magnetic susceptibility for the nearest-neighbor model for various system sizes 
and for the TPSC approach. In the latter case, the calculation is for an infinite system. The curve for 4x4 system 
is from Ref. @. 

Figure 3: Maximum in the uniform magnetic susceptibility for the next-nearest-neighbor model for 8x8 system 
and for the TPSC approach on an infinite lattice. U = results are also presented for 8x8 and for infinite size to 
estimate finite-size effects. 

Figure 4: Magnetic structure factor at the antiferromagnetic wavevector, compared with the TPSC approach in 
the infinite-size limit. The antiferromagnetic correlation length for the next-nearest-neighbor model is largest at 
half-filling but remains finite for the parameters considered here. 

Figure 5: Finite-size effects in the non-interacting uniform susceptibility when t' = —0.35. a) and b) calculations 
are at (3 = 5 for, respectively, a converged lattice and an 8 x 8 lattice, c) and d) are for the low temperature limit 
(3 = 100. 

Figure 6: Magnetic structure factor measured along the sides of the irreducible Brillouin zone illustrated in the 
insets. The TPSC calculation is for the infinite lattice, a) Half-filled case b) Filling near the van Hove singularity. 

Figure 7: Reproduced from Fig. 6 of Ref. |l7]]. Singlet d— wave corrrelation function plotted as a function of filling 
in the nearest-neighbor (t' — 0) model for U = 4, (3 = 6. The T— matrix result is the same as the U = result and is 
calculated on the same lattice sizes at the same temperature j3 — 6. a) 4 x 4 lattice, b) 8 x 8 lattice. 

Figure 8: Order parameter correlations for extended .5- wave superconductivity, a) 8 x 8 system for (n) < 1, b) 
6x6 system for the whole range of fillings. 

Figure 9: Order parameter correlations for d-wave superconductivity, a) 8 x 8 system for (n) < 1, b) 6 x 6 system 
for the whole range of fillings. 
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